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ABSTRACT 

The efficient use of Type la supernovae (SNIa) for cosmological studies requires knowledge of 
any parameter that can affect their luminosity in either systematic or statistical ways. Observa- 
tional samples of SNIa commonly use the metallicity of the host galaxy, Zh os t, as an estimator 
of the supernova progenitor metallicity, Zi a , that is one of the primary factors affecting SNIa 
magnitude. Here, we present a theoretical study of the relationship between Zi a and Zhost- We 
follow the chemical evolution of homogeneous galaxy models together with the evolution of 
the supernova rates in order to evaluate the metallicity distribution function, MDF(AZ), i.e. 
the probability that the logarithm of the metallicity of a SNIa exploding now differs in less 
than AZ from that of its host. We analyse several model galaxies aimed to represent from 
active to passive galaxies, including dwarf galaxies prone to experience supernova driven out- 
flows. We analyse as well the sensitivity of the MDF to the most uncertain ingredients of our 
approach: IMF, star-formation law, stellar lifetime, stellar yields, and SNIa delay-time distri- 
bution (DTD). Our results show a remarkable degree of agreement between the mean Zja in 
a galaxy and its Zh os t when they both are measured as the CNO abundance, especially if the 
DTD peaks at small time delays, while the average Fe abundance of host and SNIa may differ 
up to 0.4-0.6 dex in passive galaxies. The dispersion of Zia m active galaxy models is quite 
small, meaning that Zh os t is a quite good estimator of the supernova metallicity. Passive galax- 
ies present a larger dispersion, which is more pronounced in low mass galaxies. We present a 
procedure to generate random SNIa metallicities, given the host metallicity. We also discuss 
the use of different metallicity indicators: Fe vs. O, and gas-phase metallicity vs. stellar metal- 
licity. Finally, the results of the application of our formalism to a galactic catalogue ( VESPA) 
suggest that SNIa come, in average, from small metallicity progenitors both at low redshifts 
(contrary to expectations) and in galaxies with high star-formation activity. In spite of large 
uncertainties in the metallicities derived from the catalogue, the gross trends of Zi a vs. Zh 0s t 
obtained from VESPA for different galaxy types are roughly consistent with our theoretical 
estimates. 
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1 INTRODUCTION 

Metallicity is one of the few progenitor attributes that can af- 
fect the luminosity of Type la supernovae (SNIa), with impor- 
tant consequences for their use as cosmologi c al standard candle s 
jPommguez et al J 1200 ll: iTimmes et alj|2003l : [Kasen et alj|2009h . 
Recently, Bravo et al. 1 20ld) analysed the metallicity as a source 
of dispersion in the light-curve width-luminosity relationship and 
found that deriving supernova luminosities from light curve shapes 
without accounting for the supernova metallicity might lead to 
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systematic errors of up to 0.5 mag. From an observational point 
of view, the dependence of SNIa luminosities on metallicity 
might originate a systematic error of ~ 9 % in the measurement 
of the dark energy equation of state w dGallagher et al.l 120081 
ISullivan et al.ll2010f). compar able to the current statistical uncer- 
tainties. ISullivanetailfeOlOh proposed to include the host proper- 
ties (galaxy mass and metallicity) into the SNIa light-curve fitters 
in order to correct for these systematic errors. 

Up to now, attempts to measure the metallicity directly from 
superno va observations h a ve been scarce and their results un- 
certain dLentzetal.ll2000l : iTaubenberger et al.l 120081) . Measuring 
the metallicity, Z, from the X-ray emission of supernova rem- 
nants is a promising alternative but as yet has been only applied 
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to a single supernova dBadenes et al.l l2008f) . although there are 
prospects to extend this analysis to further remnants dBadenes! 
1201 d : lYamaguchi & KovamdEncl) . An alternative venue is to es- 
timate the supernova metallicity as the mean Z of its environment 
dBadenes et alj|2009h - this can be done for statistically significant 
samples of SNe in nearby galaxies by taking advant age of the long 
evolu tionary timescales of their supernova remnants dBadenes et al.l 

lioToh . 

Excluding these studies based on individual objects or super- 
nova remnants, the vast majority of our knowledge about the metal- 
licity of Type la SN progenitors has been asse mbled by measuring 
the bulk properties of their entire host galaxies. lHamuv et alj J2000h 
looked for galactic age or metal content correlations with SNIa lu- 
minosity in star-formi ng hosts, but their results were ambiguous. 
lGallagheretal] ( l2005l) studied the correlations between SNIa prop- 
erties and galaxy metallicity or, more precisely, its oxygen abun- 
dance as determine d from emission li nes in the integrated spectra of 
star-forming hosts. lEllis et all d2008l) looked for systematic trends 
of SNIa UV spectra with metallicity of the host galaxy, and found 
that the spectral variations were muc h larger than predicted by the- 
oretical models. IPrieto et all d2008l) studied the metallicity of all 
the star-forming galaxies in the Sloan Digital Sky Survey (SDSS) 
catalogue with redshifts in the range 0.01 < z < 0.04 that hosted 
SNe of some kind and noted that SNIa occur in a wi de range of 
metallicities, as low as 0.25 solar. ICooper et al l d2009l) , using data 
from the SDSS and Supernova Survey conc luded that prompt SNIa 
are more luminous in m etal-poor systems. iGallagher et al.l d2008h 
and lHowell et~ai] d2009l) . using different methodologies to estimate 
the metallicity of SNIa hosts, arrived to opposite conclusions with 
respect to the dependence of supernova luminosity on Z, although 
these discrepancies might have been solved recently dKellv et al .1 
l20ld) . We note that most of the previous work has focused on 
comparisons between the empirical observationally-determined re- 
lationship between the metall icity of hosts and t heir SNIa lumi- 
nosities on the one hand (e.g. IHowell et alj|2009h and theoretical 
predictions of correlations between the supernova meta llicity and 
luminosity on the other hand (e.g. iTimmes et al.l 120031) . thus as- 
suming implicitly that both metallicities are somehow tied. 

Estimating the metallicity of an entire galaxy is a complex 
task, involving many different observational and theoretical chal- 
lenges. On the observational side, different tracers exist for both 
gas-phase and stellar m etallicities, each with their own advan- 
tages and disadvantages l lBinnev & Merrifieldll 19981) . On the the- 
oretical side, complex modelling of the observations is often re- 
quired to arrive at a numerical estimate of the metallicity, which 
introduces biases an d uncertainties that depend on the s pecific 
technique used (e.g. IConrovetal1l2009l) . IGallagher etahl d2008h 
measured the host metallicity from absorption lines of Fe at two 
wavelengths, thus testing the gas-phase in early-type galaxies, us- 
ing single-age single-metallicity stellar population models. Then, 
to derive metallicity, they assumed s olar elementa l ratios in the 
galactic interstellar medium (ISM). IHowell et al. I d2009h deter- 
mined the host metallicity indirectly. In a first step, they deter- 
mined the host mass from spectral energy distribution fits to galaxy 
photometry. Next, they applied th e well-known mass-metallicity 
rela tionship I Tremonti et al.l |2004|) . also valid for the gas-phase. 
The lTremonti et al.l d2004l) mass-metallicity relationship was deter- 
mined in turn through measures of the nebular emission li nes of 
oxygen in star-forming galaxies. Icid Fernandes et all d2005l) com- 
pared the stellar metallicities to those of the gas-phase (measured 
through the oxygen abundance) for several thousand galaxies in 
the SDSS. They found that the relationship between both metallic- 



ities is non-linear, with a shallower than linear dependence of the 
gas-phase oxygen abundance on the stellar metallicity, and a large 
scatter of data points. 

The question we want to address in this work is: What is the 
relationship between the metallicity of a SNIa and that of its host 
galaxy? Given the complexity of the problem we want to address 
it from a statistical point of view, aiming to determine a distribu- 
tion function for the difference of metallicities between the super- 
nova and its host. Our approach is based on calculations of galac- 
tic chemical evolution together with the evolution of the super- 
nova rate. In order to understand which the best estimators of SNIa 
metallicity are we study the correlations between SNIa metallicity 
and different observational tracers of host metallicity: Fe vs. O in 
the gas-phase and stellar vs. gas-phase metallicity. In order to keep 
simplicity, we use homogeneous one-zone galactic chemical evo- 
lution models, disregarding spatial variations within each galaxy. 
Taking into account the spatial variations of metallicity would in- 
crease the scatter and corrections modelled in this paper. 

The plan of the paper is as follows. In Section|2] we describe 
the theoretical framework we adopt in order to derive the statistical 
properties of the metallicities of SNIa compared to that of their host 
galaxies. In the next section, we analyse the results, with particular 
emphasis on the reliability of different host metallicity measures 
as estimators of SNIa metallicity. Among the tracers that we con- 
sider, we rank highest the measure of SNIa metallicity through the 
gas-phase metallicity of its host. We pay special attention to the in- 
fluence of the assumed Delay Time Distribution (DTD) of SNIa. It 
turns out that the metallicity of an average SNIa is closer to that 
of its host if the D TD peaks at small time delays, as the DTD of 
iMaoz et alj feOlObh . than if the SNIa rate follows the rate of for- 
mation of white dwarfs, as proposed by IPritchet et al. I d2008l) . In 
Section|4]we explore the possibility of using current galactic cata- 
logues, such as SDSS and VESPA, to derive the metallicity distri- 
bution function of SNIa. Finally, in Section[5]we draw our conclu- 



2 THE CHEMICAL EVOLUTION OF GALAXIES AND 
THEIR SNIA 

2.1 SNIa Metallicity Distribution Function 

We define the metallicity probability density of SNIa, /"(AZ) as the 
probability that a SNIa comes from a progenitor whose metallicity 
logarithm differs from the logarithm of its host metallicity, 



AZ = log (Z Ia /Z h0St ) 



(1) 



Similarly, we define the metallicity distribution function (MDF) of 
SNIa as the cumulative metallicity probability, i.e. the probability 
that the difference between the logarithm of the metallicity of the 
SNIa progenitor and that of the host is smaller than AZ (note that 
this means AZ larger in absolute value in the usual case that the 
SNIa metallicity is smaller than the host's one), 



MDF (AZ) 



£ 



P' (AZ) dAZ . 



(2) 



We give additional details about the procedure we use to compute 
P' (AZ) at the end of Section |23l We further define the "metallicity 
correction", AZ, as the mean value of AZ, i.e.: 



AZ = log (Z Ia /Z host ) 



(3) 
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The present rate of SNIa explosions in a given host, Ri Jf), is 
the convolution over time of the Star Formation Rate (SFR), S , and 
the DTD, *F: 



Ru = f S(t' 
Jo 



) v P(r - t')dt' . 



(4) 



The probability that a SNIa that explodes now (at time f) comes 
from a progenitor born at time f = t - r is: 



PiAT) = S{t-ryV(r)IR l ,(t), 



(5) 



where r is the time spanned since progenitor birth to the moment 
the host is being observed, which we will call look-back time, and 
t is the cosmic time. Finally, the knowledge of Pu vs. time and the 
metallicity evolution, Z(f), allows to construct the MDF of SNIa. 

We have used two diffe rent prescriptions fo r the DTD: those 
of lPritchetetai] d2008l) and iMaoz etafl d2010bl in the following 



MlOb). As we will show later in Section [331 the DTD is the most 
influ ential ingredient for the determination of the MDF of SNIa. 

IPritchet et all d2008l in the following P08) found that the ob- 
servational rate of SNIa is a constant fraction, ~ 1% of the stellar 
death rate and independent of the SFR. Later, iRaskin et al.l d2009t) 
used a different method to confirm the results of P08, i.e. the rate 
of SNIa traces the rate of formation of white dwarfs (WD), with a 
constant efficiency factor: 



f S(f 



)<p{M r ) 



dM\ , 



(6) 



where <f> = (dN/dM) is the Initial Mass Function (IMF), (^) is 
the inverse of the derivative of the lifetime of a star of mass M T 
whose lifetime is precisely t, and the integral extends down to the 
evolution time of the maximum mass of a star that ends as a WD, 
M max . In this formulation, the probability that a SNIa has an age 
(i.e., its progenitor's age) r is: 



PlAT)ocS(t-T)tf>(M r ) 



dM 
If 



(7) 



On the other hand, when the DTD of Ml Ob is applied the prob- 
ability that a SNIa has an age r is obtained right from Eq.[5] 



^.(T)ocS(f-T)T(T). 



(8) 



2.2 Basic ingredients of the chemical evolution model 

We adopt a Salpeter IMF dSalpeteJ 19551) with a small-mass cut-off 
of M bo t = 0.07 M© and a high-mas s cut-off of M m = 50 M p, and 
the stellar lifetime law proposed bv lTalbot & Arnettl dl97ll) : 



t(M) = 12AT 2 ' 5 +0.005, 



(9) 



with mass in Mq and time in Gyr. In Section |2~5l we ex plore the 
sensitivity of our re sults to other choices for t he IMF (Chabrieil 
l2003l : lKrouplll2007l) and the stellar lifetime law dBuzzonill2002h . 

One of the main ingredients necessary to determine the MDF 
is the prescription for the stellar yields, includin g supernova pro- 
cessin g. Here, we adopt the yields provided by Kobav ashi et al.l 
d2006l) for stars from 13 to 40 Mq and metallicities from to 0.02. 
Stars in the above mass range are the main production site of the 
CNO elements. The lKobavashi et aid l200d) nucleosynthesis calcu- 
lations take into account the entire lifetime of the stars, from zero- 
age main sequence to the supernova or hypernova explosion, and 
provide the mass ejected in the explosion. We assume that the mass 
lost by the star during its hydrostatic evolution returns to the ISM 
the same metal fraction with which the star was born. 



The yield tables of lKobavashi et al] d2006t ) provide the ejected 
mass of each element, nij, as a function of the star mass, M,, 
and metallicity, Z,: ntj = mj(M h Zj). In most of our calculations 
we interpolate the yields as a function of metallicity, according 
to the metallicity of the ISM at the time of stellar birth. As in 
lNomotoetaI]d2006l) , we adopt a fraction of hypernovas £ hll = 0.5. 

In order to simplify, we apply the instantaneous recycling ap- 
proximation to massive stars. We integrate the yields between a 
minimum mass of 10 Mq and a maximum mass of M up = 50 Mq, 
thus obtaining a mean yield, yj, In this work, we use two measures 
of the metallicity, Z, either the CNO elements, whose mean yield is 
ycNO> or all the elements with atomic number larger than five (we 
will refer to this last group generically as 'Fe' in the following, and 
they will be our default choice for the measure of metallicity, Z, 
unless otherwise stated), whose mean yield we generically denote 
as yFiQ The integration is performed numerically by dividing the 
full mass interval in bins cen tred in the mass of each of the models 
present in the tables given bv lKobavashi et al.ld200d) . then weight- 
ing the respective yield by the IMF of mass M, and multiplying by 
the width of the mass interval. 

With respect to stars of mass below 10 Mq, we assume that 
the difference of mass between the zero-age main sequence star 
and its end-of-life compact remnant is returned to the ISM with 
the same fractional composition of metals it had at birth. Both, the 
fraction of the mass of a star that is returned to the ISM during 
its hydrostatic evolution, E\,(M), and the total fraction of the mass 
of a star that is r eturned to the I SM (e xcluding SNIa), E(M), have 
been tak en from ICatalan et all d2009h for stars of M < 10 Mq, 
and from iKobavashi et alj d200q) for stars of larger masses. Stars 
ending as SNIa return an additional mass to the ISM equal to the 
Chandrasekhar mass and composed 100% of metals, with negligi- 
ble amounts of CNO elements. 



2.3 Galactic models with mass inflow 

To study the statistical properties of the MDF of SNIa we include 
into our chemical evolution models both a time-dependent SFR 
and an infall law that account for the interaction with the galac- 
tic halo. For the SFR w e adopt the Schmidt law as implemented by 
IKobavashi et alj|200d) : 



S(f) = mM)/T sfr , 



(10) 



where T str is the timescale of star formation, and m g (r) is the mass 
of gas in the ISM at time t. In the litera ture, ther e are o ther pre- 
scriptions for t he SFR law. For instance. ISandagel (fl986l see also 
iGavazzi et ail d2002l) : iGreggio & Cappellarol d2009l) ) proposed to 
use a universal exponential SFR as a function of time: 



Sit) oc — exp -(r/T sfr ) 2 



(11) 



with r s f r ~ 2 - 4 Gyr for early-type galaxies, r s f r ~ 10-13 Gyr 
for spirals, and T s f, ~ 25 Gyr for dwarf irregulars. However, the use 

1 For SNIa theory both measures of metallicity are meaningful. On one 
side, the total neutron excess in the exploding white dwarf, related to the 
overall metallicity of the progenitor, plays a role in the timescale of electron 
captures that control the fraction of incinerated mass that goes to radioactive 
56 Ni and empowers the light curve. On the other side, the CNO elements 
present in the progenitor turn into 22 Ne in the white dwarf, which has a rel- 
evant role as a catalyst of the thermonuclear reactions in SNIa a nd may also 
influe nce the brightness of the supernova (for a discussion see lBravo et alj 
lioiol and references therein). 
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of the Sandage SFR law can lead to non-positive gas mass if the 
constant of proportionality is not accurately chosen and/or other 
ingredients such as inflow/outflow are present in the model. It is 
for this reason that we prefer to adhere to the Schmidt law of SFR. 

To include the infall of m aterial from the galactic halo we fol- 
low the conservative model of iKobavashi et al.l and assume 
that the total mass of the system formed by gas, stars and halo is 
constant and equal to Mw Initially, the whole mass is in the halo, 
and the mass infall rate as a function of time is: 
M b0 



A i 



pat. 



(17) 



F(f): 



■ exp [-f/T mf ] 



(12) 



where T inf is the timescale of infall. IKobavashi etafl feOOCl) propose 
values of r- mf = 0.1 Gyr as representative of elliptical galaxies, 
10.9 Gyr to account for Galaxies of type Sbc/Sc, and 60.5 Gyr to 
represent types Scd/Sd. In this work, we use T sfr and r inf as free 
parameters in order to obtain models for different types of galaxies 
The set of differential equations to integrate in these models is 
the following one (in addition to Eqs.[7]or[8): 



-S(t) + F(f)+R Ia M, 



Ch 



(f>(M)S (t - r M )E{M)MAM , 

M(t) 



-S(t)[Xj{f) -y,]+ R la M ch 

+ (p(M)S(t-T M )E h (M)Xj(t-T M )MdM , 

Jm(i) 



(13) 



(14) 



where M(t) is the inverse of the stellar lifetime function (Eq. |9), 
Xj = nij/m g is the mass fraction of elements j in the gas at time f, 
tm is given by Eq. [9] 5 a function of time through m E (t), the sub- 
script j in Eq. [14] stands for either 'CNO' or 'Fe', and i?i a is com- 
puted using Eq.|6](for the P08 DTD) or Eq.@](for the MlOb DTD). 
In this equation, we assume that the infall gas is of primordial com- 
position and does not contribute to the mass of metals. Seemingly, 
the term R^Mq^ is only added to the equation of evolution of m Fc , 
but not to that of m mo . The last term in Eq. [14] accounts for the 
fraction of metals in the protostellar nebula that returns to the ISM 
during the hydrostatic stellar evolution. 

If the mass fraction Xj were a monotonic function of time, 
the metallicity probability density might be obtained directly by 
combining Eqs.[5][T3] and !14FI 



dAZ 1 Xj- 



(15) 



where we have used the definition of AZ (Eq.QJ and the relationship 
between and Ri a , Pi a = d^i a /df, Xj is the measure of Z, and 



nij m.j 

X: = -m g . 

m„ mi 



(16) 



However, it is not guaranteed that the metallicity is a monotonic 
function of time, so we have to resort to a different procedure to ob- 
tain the metallicity probability density. In this work we have com- 
puted P' (AZ) through a Monte Carlo calculation. First, we solve 
Eqs.[5][T3] and [14] to obtain P Ia and Xj as functions of time. Sec- 
ond, we generate randomly the distribution of birth times of SNIa 
according to Pi a , by drawing N random numbers, xu from a uni- 
form distribution and assigning them birth times t t such that 

2 A normalization constant is missing in Eq. [15] We have omitted it to 
improve readability. 



For each t, thus obtained, the corresponding AZ is identified and 
the probability density is easily calculated. For this last calculation 
we have used 100 bins in AZ. The results shown in this paper were 
generated with N = 100000 random numbers per each j each time 
the chemical evolution of a galaxy was computed. 



2.4 Models of passive galaxies with mass outflow 

Elliptical galaxies span large luminosity and mass ranges, their 
light emission is domin ated by red giant stars, and contain no gas 
(for a recent review see lMatteuccill2008[ and references therein). 
Ellipticals are metal rich, although t heir mean stellar metallicity i s 
in the ranee [Fe] = from -0.8 to +0.3 jKobavashi & Arimotol 19991) , 
in general the metallicity grows with the mass of the galaxy. The 
history of star formation in elliptical galaxies is still controversial, 
the two main scenarios being the single starburst model (mono- 
lithic scenario) and the hierarchical model. In the former, ellipti- 
cal galaxies are assumed to have formed at high redshift through a 
short and intense starburst as a result of dissipative collapse of pro- 
toga lactic gas clouds followed by a supernova driven galactic wind 



(e.g.lLarsonll974lArimoto & Yoshiill987l : lMatteucci & Tornambd 
1 19871 ; IPipino & Matteuccil |2004|) . In the hierarchical model, el- 
liptical galaxies form at relatively recent epochs through merg- 
ers of gaseous galaxies, with continuous star formation through a 
wide redshift range (e.g. [White & Reeslll978l: baugh et al.|ll998l : 
iKauffmann & Charloti ri998; Stei nmetz & Navarroll2002ir In every 
case, elliptical galaxies should be characterized by short and in- 
tense starburst(s) during which a large number of massive stars 
form and explode in a short time interval, releasing a large energy 
into the gas leading to massive outflows that disperse the super- 
nova yields into the intergalactic medium, especially for small mass 
galaxies. 

The question of which model best reproduces the observa- 
tional constraints is not settled yet. Since we do not pretend to 
address all th e complexities of the form ation of elliptical galax- 
ies (see, e.g., IPipino & Matt eucci 200^), we have just chosen to 
adopt a simple prescription intended to reproduce the chemical evo- 
lut ion in the monol ithic scenario: the so-called "leaky-box model" 
of lHartwick! jl976l) . In this prescription, the effect of the super- 
nova driven outflo w is encapsulated in the use of an effective yield, 
>';.cff = >'i7(l+ c out) dBinnev & Merrifieldll998l) . where y, is the usual 
yield of element i, y,_ e (f is a reduced yield that goes to the ISM and 
accounts for all the material lost to outflow motions, and c out > 
is a constant. The leaky-box model allows us to simulate ellipti- 
cal galaxies of different masses, where large values of c out repre- 
sent low-mass low-metallicity early-type galaxy models. We have 
computed models of early-type galaxies with infall timescale and 
SFR timescale in the range 0.1-0.5 Gyr (see IPipino & Matteuccil 
120041 for an study of the correlation between both timescales and 
the size and luminosity of the galaxy), and have explored a range 
of Cout from (no outflow, a model for giant elliptical galaxies) to 
10 (meant to represent dwarf ellipticals). 



2.5 Sensitivity to the IMF, stellar lifetime, and stellar yields 

Our objective is to understand the relationship between the metal- 
licities of SNIa and those of their hosts. The sensitivity of this rela- 
tionship to uncertain ingredients of the theoretical model depends 
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Figure 1. Comparison of th e IMFs u s ed in the pres ent work in a ra nge of 
masses of SNIa progenitor s: ISabeterl fl95l (solid), iKroupal d2007l) (short 
dash), and|Chabriei (2003) (long dash). The vertical axis gives each IMF as 
a function of star mass normalized to the same IMF law for M = 2 Mq. 



well those of lKroupal d2007h and lChabrierl d2003l) . Figure Q] shows 
the behaviour of these IMF within the range of masses of SNIa 
progen itors . In this mass r ange, the normalized IMF of ISalpeterl 
andlChabrierl d2003l) are nearly indistinguishable, while the 
IMF o f iKroup al 120071) gives slightly less weight to the more mas- 
sive SNIa progenitors, consequently favouring slightly longer SNIa 
delay times. These small differences translate in small variations in 
the MDF of SNIa. The difference in the metallicity correction, AZ, 
between all these IMF is smaller than ~ 0.08 dex for the whole 
sample of galactic models we have generated. 

To address the dependence of our results on the stellar lifetime 
law, we have r epeated the calculations using the prescriptions of 
iBuzzonil d2002h instead of iTalbot & Arnettl fl97lh . In general, the 



Buzzonil i 200^ ) are sli ghtly shorter (- 



lifetimes obtained followin 

50%) than those given by ITalbot & Arnettl Jl97lT) . Although this 
leads SNIa to follow tighter the SFR (as we are here working with 
the DTD of P08), the metallicity correction does not change by 
more than ~ 0.07 dex between both lifetime laws. 

In order to test the sensitivity of our results to the prescrip- 
tion of the stellar yi elds we have repeated th e calculations with the 
yields ntj given bv lKobavashi eUri1 f200ll) fixed at two extreme 
metallicities: either Z, = or Z, = 0.02. The result is that the MDF 
of SNIa is insensitive to the choice of stellar yields. This insensitiv- 
ity is a consequence of our definition of the MDF as a function of 
the ratio of Zi a to Zh ost . Hence, even if the absolute values of both 
metallicites change with the stellar yield prescription, their ratio is 
extremely insensitive to them. 



on the way these ingredients affect the host chemistry in compar- 
ison to the SNIa rate. Here we use a simplified treatment of the 
galactic chemical evolution that allows obtaining MDF indepen- 
dent of the SFR, and the DTD of P08. 

We start by considering the simultaneous evolution of metal- 
licity and SNIa rates in a closed-box model. We further simplify the 
model by taking a constant SFR, So, through the complete galac- 
tic life, and assuming the instantaneous recycling approximation 
for stars more massive than 10 Mq. Within these assumptions, the 
evolution of the mass of gas and of the mass of CNO elements in 
the gas phase, «i C no, are given by the following differential equa- 
tions: 



Jmu) 



</>(M)E(M)MdM 



- -So [Xcno(0 - VCNo] , 



(18) 
(19) 



where vcno incorporates now the last term in Eq. [14] Within the 
same assumptions, Eqs.|7]and[8]simplify to: 



(20) 



(21) 



and the MDF can be obtained formally from: 

dtfla Y(T) 

oc , 

dmcNO XcnoW - ycNO 

using the results of the integration of Eqs. [T8]and[T9] The last ex- 
pression shows that the MDF thu s obtained is ind ependent of So- 

The IMF we use is that of lSalpeterl dl955l) . The results we 
obtain are insensitive to the precise values of the cut-off within the 
range we have explored: M bat = 0.07 - 0.1 Mq and M up = 50 - 
80 Mq. 

Besides the IMF of lSalpeteil dl955T) . we have considered as 



3 THEORETICAL RESULTS 

In this section, we will present the results obtained with the models 
including inflow and the models of passive galaxies. Unless other- 
wise stated, we solve Eqs.[T3]and[T4]from t = to / = / ga | = 13 Gyr. 



3.1 Prototype galaxies 

We start by analysing the behaviour of typical galaxies meant to 
represent from ellipticals to late-type galaxies. Table Q~]gives their 
main parameters together with the final mass fractions of gas and 
stars with respect to the total mass of the system, the final metallic- 
ity (both as given by all metals and by CNO elements) of the host, 
and the statistical properties of the corresponding MDF of SNIa: 
metallicity correction and standard deviation of the distribution, a. 
We computed all the model galaxies discussed in this subsection 
with no outflow. Mo dels A to C are i n goo d agreement with the 
models in Table 2 of iKobavashi et all d200d ) for galaxies of types 
E, Sc, and Sd (the remaining mass of gas and the mass of stars at 
t = ; ga i of the models, as well as the final metallicities agree reason- 
ably well given that the model of evolution of the rate of SNIa we 
use are not the same as their). The SFR of model D is more bursty 
(Fig. 0, and we include it here as a representative of the early- 
type galaxy models used in the rest of the paper (even though it has 
c ou t = 0, which would correspond to a massive early- type galaxy). 
In this and the next subsection, we present the results obtained with 
the DTD of P08, while we delay to Section [331 the discussion of 
the differences with respect to the DTD of Ml 0b. 

In Fig. [3] we show the distribution of birth-times for a SNIa 
exploding at / = t m \. The curves peak at the present time because 
the DTD is strongly biased towards small delay times. This ef- 
fect is more pronounced in late-type galaxies, for which the star- 
formation activity is either slightly decreasing or increasing. Early- 
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Table 1. Parameters of some typical galaxy models. 



Pritchet-DTD models Maoz-DTD models 





A 


B 


c 


D 


A' 


B' 


C 


D' 


Ttaf (Gyr)" 


0.1 


9.0 


50.0 


0.1 


0.1 


9.0 


50.0 


0.1 


r sfr (Gyr)* 


4.0 


2.0 


2.0 


0.1 


4.0 


2.0 


2.0 


0.1 


DTD^ 


P08 


P08 


P08 


P08 


Ml 0b 


MlOb 


MlOb 


MlOb 


f d 


0.09 


0.09 


0.04 





0.09 


0.09 


0.04 







0.91 


0.67 


0.19 


1 


0.91 


0.67 


0.19 


1 




0.029 


0.015 


0.012 


0.036 


0.041 


0.022 


0.017 


0.035 


AZ (dex£ 


-0.53 


-0.15 


-0.089 


-0.54 


-0.21 


-0.047 


-0.029 


-0.33 




0.49 


0.19 


0.14 


0.35 


0.39 


0.11 


0.071 


0.41 


X(CNO) h0St i 


0.017 


0.0091 


0.0076 


0.0085 


0.017 


0.0091 


0.0076 


0.0085 


AX (CNO) (dex}£ 


-0.48 


-0.12 


-0.070 


-0.13 


-0.19 


-0.024 


-0.017 


-0.093 


cr{ AX (CNO) £ 


0.47 


0.18 


0.13 


0.31 


0.36 


0.10 


0.062 


0.28 



a Timescale of infall 

^ Tim escale of star formati on 

c P08: |Pritchet et al.1 120081) ; M10b: lMaoz et alj (2010bl) . 

d Fraction of gas mass in the ISM at t = 13 Gyr with respect to the initial baryonic halo mass, M;,o. 
e Fraction of stellar mass at ( = 13 Gyr with respect to the initial baryonic halo mass, M/,o. 
f Gas metallicity (mass fraction of all metals) at t = 13 Gyr. 
8 AZ^toiCZia/Zhost). 

^ Standard deviation of the distribution of SNIa metallicities, MDF. 
' Mass fraction of CNO elements in the gas at t = 13 Gyr. 
J AX (CNO) = ioi [X (CNO) Ia IX (CNO) host ] . 

* Standard deviation of the distribution of CNO mass fractions in SNIa. 



type galaxies display a similar peak at quite early times, belong- 
ing to the epoch of intense star-formation activity, but the impor- 
tance of this initial peak varies according to the duration of the 
star-formation epoch. Thus, SNIa from early-type galaxy models 
come from two different populations whose properties can be quite 
different: an old, predominantly low-mass, population belonging to 
the star-formation peak, and a young population characterized by 
larger masses. 

Figure [4] shows the MDF of SNIa, for the same typical galax- 
ies as those in the left panel of Fig. [3] as a function of AZ. As 
can be deduced from Fig. [3] the younger the stellar population of a 
galaxy is (as in models B and C) the smaller the scatter in metallic- 
ity of potential SNIa progenitors. The late-type galaxy distribution 
is strongly biased towards a zero correction and display a short tail, 
down to AZ - — 0.2 dex for a galaxy that is still in its initial phases 
of star formation (T illf = 50 Gyr, meant to represent an irregular 
galaxy), and AZ ~ -0.4 dex for a galaxy that has already gone 
through most of the infall process (r inf = 9 Gyr). The probability of 
a SNIa metallicity differing in less than 0. 1 dex (in absolute value) 
from Zh 0S t is 75% for the irregular galaxy model and 57% for the 
model with T inf = 9 Gyr. 

Early-type galaxies show a broader MDF implying that SNIa 
from these hosts will have a larger intrinsic metallicity dispersion, 
even if they come from the same host! Again, this behaviour can be 
understood from the probability density of SNIa progenitors shown 
in Fig. [3] Models A and D have a substantial contribution of old 
potential SNIa progenitors, thus they show a large scatter in their 
MDF. For the model with r illf = 0. 1 Gyr and T sfr = 4 Gyr half of the 
SNIa have a metallicity differing more than 0.38 dex (in absolute 
value) from Zhost. while for the model with Tin! = 0. 1 Gyr and T s f, = 
0. 1 Gyr the median of the correction is -0.49 dex. 

Figure[5]shows the MDF of SNIa when the metallicity is mea- 



sured as the abundance of CNO elements. The main difference with 
the MDF derived for 'Fe' is found in the behaviour of the older 
galaxy model, i.e. the one with T- m { = 0. 1 Gyr and r s f, = 0. 1 Gyr. 
The reason of this behaviour can be understood by looking at Fig. [3] 
(bottom right panel), where we display the evolution with time of 
the mass fraction of CNO elements in the gas. The long-dashed 
(magenta) curve shows the evolution following the short burst of 
star formation coincident with the infall epoch. Initially the abun- 
dance of CNO elements grows fast due to the self-enrichment of the 
gas due to efficient reprocessing driven by massive stars. Once the 
star formation burst is over, there are no more massive stars present, 
and intermediate and low mass stars enrich the interstellar medium 
at the end of their lifetime with yields representative of their birth 
time, on average less rich on CNO elements than at the peak of the 
curve. This causes the CNO abundance to decrease rapidly until it 
reaches a stable value at about 60% of the peak CNO mass frac- 
tion. Therefore, the SNIa that explode at / = 13 Gyr can have birth 
metallicites larger than that of the host, as seen in Fig. [5] 

Because SNIa are observed at different redshifts, it is also in- 
teresting to know if their MDF varies with the elapsed time of 
galactic evolution, f gal . We have integrated again the evolutionary 
equations with f eiu = 11 and 13.7 Gyr. Within this range of t g a the 
"metallicity correction" changes by less than ~ 10%. 



3.2 Statistical properties of the SNIa MDF for different types 
of galaxy 

In this section, we use a sample of galactic models with varying r sfr , 
Tinf, and c out to study the statistical properties of the relationship 
between Z Ia and Z host . We have built models for early-type galaxies 
by using small values of the timescales associated to infall and SFR: 
either r- m { = 0. 1 Gyr or 0.5 Gyr, and either r s f r =0.1 Gyr or 0.5 Gyr, 
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Figure 3. Chemical evolution of the prototypical galaxies (the same models as in Fig.[2]and Table[TJ. Top row: Probability density that a SNIa exploding at 
r = 13 Gyr was born at a given cosmic time obtained with the DTD of P08 (left), and with the DTD of MlOb (right). The infall of gas that primes the formation 
of stars begins at a cosmic time t = for all the models. Bottom row: Evolution with time of the gas-phase metallicity (left), and the mass fraction of CNO 
elements (right). 



and with different degrees of outflow, controlled by the parameter 
c out . Thus, for each value of c out we show four models, obtained 
with the aforementioned combinations of r m f and T s f r . In general, 
the location of these models in the plane Z Ia vs Z hosl is driven by 
Cout, while the variation with tw and r sfr is small. 

In order to generate models for late-type galaxies-^ we left the 
infall and SFR timescales vary in the ranges: Tmt = 0. 1 : 60 Gyr 



2 In the following we will refer as late-type galaxy models those models 



and T sfr = 1:6 Gyr. In order to avoid too extreme or unrealis- 
tic galaxy models we have further constrained our sample to have 
0.05 < /gaJ/sta < 3 and -0.6 < log (Z hosl /Z ) < 0.5. 

Figure [6] (left panel: DTD of P08) shows the metallicity cor- 
rection, to be applied to Z host in order to estimate Z id , as a function 
of the host metallicity. The two kind of galaxies show different be- 



generated with r s f r > 1 Gyr and c out = 0. Note that this includes models A 
to C and A' to C in TableQ] 
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Figure 2. Star formation rates as a function of time for four galaxy mod- 
els in Table[T] The different curves are meant to represent typical galaxies: 
model D in Table[T](long-dashed magenta) belongs to an early-type galaxy, 
while models B (dotted black) and C (short-dashed blue) are representative 
of late-type galaxies, and model A (solid red) is a transition galaxy char- 
acterized by a short infall timescale and a moderately long star formation 
rate timescale. The vertical axis gives the SFR normalized by the total mass 
of the system, Mho (Mho is also the mass of the halo at t = and should 
be much smaller for irregulars and dwarf ellipticals than for normal-size 
galaxies). 
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Figure 4. Metallicity distribution function (MDF) of SNIa as a function of 
the metallicity correction with respect to Z(, ost , for the same models as in 
the left panel of Fig. [3] 



Figure 5. Metallicity distribution function (MDF) of SNIa as a function 
of the metallicity correction with respect to the host metallicity, but taking 
the CNO abundance as a measure of metallicity. The models shown are the 
same as in the left panel of Fig. [5] 

haviour: while there is an approximately linear dependence of the 
correction on Z hosl for late-type galaxies, the correction is about 
constant for early-type galaxies of every metallicity, i.e. the disper- 
sion of the corrections is much smaller than the range of host metal- 
licities. The linear function that fits the late-type galaxies passes 
as well through the early-type ones with no outflow (c out = 0). As 
c oul increases, the final metallicities of the early-type hosts decrease 
while the metallicity correction remains the same. 

Figure [7] shows the metallicity correction belonging to CNO 
abundances. It displays the same gross properties as those seen in 
Fig. [6] with one important difference concerning early-type galax- 
ies: their metallicity corrections are now among the lowest of all 
models displayed. The reason for this behaviour is again the pecu- 
liar MDF of CNO elements in SNIa from this kind of galaxies as 
seen in Fig. [5] The dispersion of their MDF is, however, similar to 
that of galaxies with a much larger metallicity correction (compare 
models A and D in Table[T]l. 

In many observational studies, it is necessary to generate a ran- 
dom set of SNIa with stochastic metallicities in order to compare 
their statistical properties with those of the observational sample. 
The linear fit of AZ vs. [Zh 0s t] shown in Fig. [6] together with the 
MDF may be used to generate such a random set of SNIa, with 
stochastic metallicities related (but not equal) to their host metal- 
licities. The simpler way to do this is to approximate the MDF of 
SNIa as a straight line in the semi-logarithmic plot shown in Fig.|4] 

log(MDF) = o-AZ , (22) 

that determines the probability density, 

/"(AZ) = a exp (aAZ) , (23) 
the mean metallicity correction, 

AZ = -l/a, (24) 
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Figure 6. The metallicity correction, i.e. the mean difference between the logarithms of Zi a and the host metallicity, Z(, ost , as a function of the last. Blue circles 
represent early-type galaxies, which we have built with small values of the timescales associated to infall and SFR: Tj n f = 0.1 and 0.5 Gyr, and r s f r = 0.1 and 
0.5 Gyr, and with different degrees of outflow, controlled by the parameter c ollt ■ In this and coming figures we have plotted as late-type galaxies (red crosses) all 
the models with r s f, > 1 Gyr. The error bar at the top left of each panel is the typical uncertainty in observationally determined [Zh ost ] (see Section l3~4t . Left: 
Models computed with the DTD of P08. The solid red line is a linear fit to the data points belonging to late-type galaxies given by: AZ = —1.22 [Zh 0s t] ~ 0.28, 
where we have adopted a solar metal fraction Zq = 0.01895. The dashed lines are linear fits belonging to early-type galaxies (blue) and late-type galaxies 
(red) in the VESPA catalogue computed using the DTD of P08 (see Section |4~2V The filled symbols belong to models C, B, A, and D in Table[T](from top 
to bottom). Right: Models computed with the DTD of MlOb (note the different range of the vertical axis as compared to the left panel). The solid red line is 
a linear fit to the data points belonging to late-type galaxies (red crosses), given by: AZ = -0.54 [Z(, ost ] - 0.03. The dashed lines are linear fits belonging to 
early-type galaxies (blue) and late-type galaxies (red) in the VESPA catalogue computed using the DTD of MlOb (see Section l4~2V The filled symbols belong 
to models C\ B', A', and D' in Table[T](from top to bottom). 



and the standard deviation, cr, of the distribution of AZ: 



a = l/a = -AZ. 



(25) 



In TableQ]we give the standard deviation of the MDF of our typical 
galaxy models. As can be seen, the relationship between AZ and 
cr derived in Eq. [25] approximately holds for all the models, with 
the exception of the distribution of CNO in models D and D', for 
which a simple proportionality law as that of Eq.|22]is not a good 
fit (see Fig. [5}. 

To generate a stochastic distribution of SNIa metallicites in a 
late-type galaxy of known Z host , one would have first to apply the 
linear function of Fig. [6] (or that in Fig. |7]if what is desired is the 
distribution of mass fractions of the CNO elements) to obtain the 
metallicity correction AZ and cr (= -AZ). For an early-type host 
the procedure is even simpler, as the metallicity correction takes on 
a nearly constant value for this kind of galaxies. 



3.3 Choice of the Delay Time Distribution 

The Delay Time Distribution is a property specific of SNIa that 
has rather small imprint on the galactic chemical evolution of 
CNO elements but affects strongly the temporal distribution of 
SNIa, hence it is expected to influence appre ciably the MDF. The 
DTD of SNIa is a matter of current debate ( Strolger et all 120041 ; 
Scannapieco & Bildstedl2005l lGreggidl2005l ; |Greggio et al1l2008l ; 
Maozetalj2010d) . P08 proposed a model in which the rate of SNIa 



was proportional to the WD formation rate. This model reproduced 
satisfactorily the trend of SNIa rates as a function of the specific 
SFR in galaxies of ages between 1 and 13 Gyr. In a recent work, 
iMaoz et all ( l2010al) presented a method to recover the supernova 
DTD that simultaneously takes into account supernova data and the 
reconstructed star formation history of the individual galaxies in a 
survey, and applied this method to the events present in the LOSS 
SN su rvey and in the VESPA database. Independently. iBrandt et alj 
feoich analyzed light curves and host galaxy spectra of 101 SNIa 
from the SDSS u sing VESPA algorith ms, from which they derived 
the DTD of SNIa. llvIaoz et alj l l2010d) found that their data required 
a stronger short delay-time component than allowed by the model 
of P08, i.e. they called for a higher contribution from large mass 
progenitors. The VESPA database is organized in temporal bins, 
whose resolution is better at the present epoch and degrades with 
lo ok-back time (for fu rther details see Section|4), hence the results 
of lMaoz etailfcoiod) should be more sensitive to the short delay 
time component of the SNIa DT D. Thus, the discrep ancy between 
the models proposed by P08 and M aoz et al. (2010a) might be due 
to the different sensitivities to the time ranges of the SNIa DTD 
(but see MlOb). Anyway, it is necessary to check the imprint that 
the different DTD have on the MDF of SNIa. 

Having used a model based on P08 in the previous sections, 
we now discuss the results obtained with the DTD of MlOb, Eq. 13 
in their paper, 



¥(/) = 0.7 x 10~V SN M Q ' Gyr 



(26) 
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Figure 7. The metallicity correction with respect to the gas-phase metallicity of the host, when the metallicity is measured as the CNO abundance. The 
error bar at the bottom left of each panel is the typical uncertainty in observationally determined [X(CNO)h osl ] (see Section |3~4V Left: Models computed 
with the DTD of P08 (note the different range of the vertical axis as compared to the left panel). The straight line is a linear fit to the data points belonging 
to late-type galaxies (red crosses), given by: log (X(CNO)i il /X(CNO)h st) = -1.24 [X(CNO)h os t] - 0.35, where we have adopted a solar CNO mass fraction 
X(CNO)q = 0.0134. From top to bottom, the filled circles belong to models C, B, D, and A in TablefJ] Right: Models computed with the DTD of MIOb. The 
straight line is a linear fit to the data points belonging to late-type galaxies (red crosses), given by: log (X(CNO)i a /X(CNO)host) = -0.53 [X(CNO)h os t] _ 0.19. 
From top to bottom, the filled circles belong to models C, B', D', and A' in Table[T] 



where s = -1.2, t is measured in Gyr, and we have applied a factor 
0.7 to the DTD of M 10b to convert back from their 'diet-Salpeter' 
IMF to ours. Models A' to D' in TableQ]were computed using their 
DTD in Eqs.[4]and[5] The right panel of Fig.[3]shows the distribu- 
tion of birth-times for a SNIa exploding at t = f gal for these model 
galaxies. With respect to the DTD of P08 (left panel of Fig. [3} the 
new DTD increases the contribution from young stars (small delay- 
time) with respect to old stars. Late-type galaxies are scarcely af- 
fected at all, but early-type galaxies display a quite different trend 
in which the contribution of the population belonging to the initial 
star-formation peak is dramatically reduced. With this new DTD, 
the differences between galaxy-types strongly smooth out. 

The metallicity correction obtained with the DTD of MIOb is 
shown in the right panel of Figs. [6] and [7J for both late-type and 
early-type galaxies. Overall, the correction is significantly smaller 
than that obtained using the DTD of P08, as expected given the 
larger contribution from prompt SNIa implied by the DTD of 
MIOb. 

From these results it is clear that the DTD of SNIa is the 
most important factor determining the magnitude of the correction 
needed to obtain the metallicity of SNIa from those of their host 
galaxies. We stress that any of the above DTD models (P08 vs. 
MIOb) is subject to a large degree of observational uncertainty. For 
instance, MIOb discuss two possible approaches to determine the 
power law exponent in Eq.[26] based on what they call the 'optimal- 
iron constraint' and the 'minimal-iron constraint', which lead to 
different exponents, s, in the range from -1.5 to -0.9 (in this study, 
we have adopted in Eq. [26] the mean, i.e. -1.2). It is interesting to 
consider how much would change the plots in Figs. [6] and [7] if * 
covered the full range given by MIOb. We have repeated the calcu- 



lations with s = -1.5 and s = -0.9. For the first case, belonging to 
an extreme case of the 'optimal-iron constraint', the straight-line fit 
to our late-type galaxies changes to AZ = -0.28 [Zh ost ] + 0.06 while 
the correction for early-type galaxies changes by ~ +0.12 dex. 
For the last case, belonging to an extreme case of the 'minimal- 
iron constraint', the fit of the late-type host corrections changes to 
AZ = -0.89 [Z hosl ] - 0.16, and the corrections for early-type hosts 
changes by ~ -0.04 dex. We can see that the uncertainty in the 
DTD leads to uncertainties in AZ of the same order as the own cor- 
rection. We conclude that an accurate determination of the DTD of 
SNIa is a prerequisite to achieve a good understanding of the dif- 
ferences between the metallicities of these supernovae and those of 
their host galaxies. 

3.4 Is it the host metallicity a good estimator of the SNIa 
metallicity? 

Here we discuss the reliability of different host metallicity estima- 
tors as for establishing the statistical properties of SNIa: gas-phase 
vs. stellar metallicity, Z stal , and gaseous CNO abundance vs. F^l 
abundance. 

We have evaluated the mean host stellar metallicity as the 
mean of the metallicity of the ISM at star birth, weighted by the 
SFR and by the IMF, for all the stars whose mean-sequence lifetime 
is larger than the elapsed time since their formation. The mean stel- 
lar metallicity derived by observational studies is weighted by the 
contribution of stars in different phases of their life, with weights 

3 Recall that what we are calling here 'Fe' represents all the metals. 
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depending on the particular spectral feature used for measuring the 
metallicity. Even though both ways to define a mean stellar metal- 
licity may give different values, we expect that both will follow the 
same qualitative trends. 

Fig. [8] shows the metallicity correction to apply if the host 
galaxy metallicity is measured from the stellar population (and the 
DTD of Ml Ob is used). In comparison with the trend shown in the 
right panel of Fig. [6] we notice that the corrections to apply to the 
metallicities of early-type galaxies continues being independent of 
host metallicity, but now they cluster around AZ = 0. Late-type 
galaxies display also a similar behaviour as in Fig. [6] with changed 
sign, but they could not be fitted by a linear function and show 
larger dispersion for given Z host . Thus, if possible, it is advisable to 
measure the metallicity of SNIa in late-type hosts as that of their 
gas. These results can be understood by computing the distribution 
of birth-times of present-day (at t = / sal = 13 Gyr) stars. We have 
calculated the distribution of birth times as the fraction of the mass 
in present-day stars represented by the stars born at time t that are 
still alive, 



dra s . 



dt 



J~M(f gal -() 
</>(M)MdM . 



(27) 



The distribution of birth-times follows closely the SFR of each 
model (see Fig.0 because 99% of the stellar mass at birth belongs 
to stars smaller than ~ 0.5 M Q (for our Salpeter IMF), hence the 
integral in Eq. [27] is always close to ~ 1. Comparing Fig. [2] with 
the probability density of birth times of SNIa shown in Fig. [3] (top 
right) it is easy to see that the average age of SNIa progenitors is 
much less than the average age of the stellar population, especially 
because the peak of P Ia at times close to f ga i. In galaxies charac- 
terized by a monotonic increase of the metallicity (models A', B', 
and C) the immediate consequence is that the metallicity of SNIa 
is larger than the average stellar metallicity. Early-type galaxies, 
like model D', display a non-monotonic evolution of the metallic- 
ity, that peaks later than the SFR. For instance, for model D' the 
peak of the SFR occurs at t — 0.1 Gyr, while the peak of Z gas is 
found at ~ 0.85 Gyr and that of Xcno takes place at ~ 0.60 Gyr. 
Thus, the stellar population is chiefly drawn from low-metallicity 
gas, while SNIa have an appreciable contribution from times close 
to present, at which the gas metallicity is much larger. 

In Fig. [7j we show the metallicity correction when the gas- 
phase metallicity is estimated from the mass fraction of CNO el- 
ements in the gas. In comparison with Fig. [6] the CNO metallic- 
ity correction is nearly identical to the Fe correction for late-type 
galaxies but significantly smaller for early-type hosts, for a given 
DTD. As can be seen in Fig. [3] (bottom row) the chemical evolu- 
tions of Fe and CNO elements in the galaxy models A', B', and C 
are in close agreement, thus the MDF of both element groups are 
similar (compare Fig.|4]with Fig. [5}. Thus, the use of Fe instead of 
CNO abundances does not seem to represent an improvement in the 
estimates of SNIa metallicities obtained from those of their hosts. 
The host CNO abundance of early-type galaxies is in much better 
agreement with the mean SNIa one than are the corresponding to- 
tal metal abundances, although their dispersions are similar. Using 
again the chemical evolution computed for the prototypical galaxy 
model D', displayed in the bottom row of Fig. [3] and the proba- 
bility density of SNIa shown in Fig. [3] the reason for this closer 
agreement between CNO abundances in SNIa and the gas, as com- 
pared with Fe abundances, is that SNIa from early-type galaxies 
have non-negligible contributions from quite early times and near 
present time. The CNO abundances are low at both times, because 
presently there are no high-mass stars contributing to their synthe- 
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Figure 8. The metallicity correction, taking as a reference the mean stellar 
metallicity, Z star- , instead of the gas metallicity, Zj, ost . Here we have used the 
DTD of MIOb. The filled symbols belong to models A', B', C, and D' in 
Table[T](from top to bottom). The error bar at the top left of the figure is the 
typical uncertainty in observationally determined [Z slal ] (see discussion at 
the end of Section 1X41 . 



sis. On the other hand, the Fe abundance increases continuously 
during the last 11 Gyr due to the contribution of SNIa, with the 
result that young SNIa progenitors are drawn from gas with high 
Fe abundance while old SNIa progenitors were born from gas with 
quite small Fe abundance. Hence, the average discrepancy between 
present-day gas Fe abundance and that of SNIa is larger than in the 
case of the CNO abundances. 

A meaningful evaluation of the corrections we propose to ob- 
tain SNIa metallicities from the different metallicity measures of 
their hosts can only be done by considering the observational errors 
in the galactic metallicity and in the DTD of SNIa. As discussed 
in the Introduction, galactic metallicities are measured using dif- 
ferent methods, all of th em subject to a suite of experimental and 
theoretical uncertainties. Galla gher et al.l d2008l) measured the Fe 
abundance of the gas of early-type hosts of SNIa with quite dif- 
ferent errors, depending on the host. An average error in [Z host ] of 
0.1-0.2 dex can be deduced from their Table 4 and fr om the com- 
pariso n the y perform with meta llicities measured by iTrager et al.l 
hoOOl) and lThomas et alj d2005l) for a subset of galaxies the y have 



in common. On the other hand, the iTremonti et ah |j2004 mass- 
metallicity relationship, used by iHowell et alj d2009l) to measure 
host metallicity, is sensitive to [O/H] in the gas phase. The scat- 
ter in this relationship is of order ~ 0.1 dex. Finally, the obser- 
vational determinati on of [Z star ] is subject to s ystematic biases of 
order 0.1 - 0.2 dex dCid Fernandes et al.ll2005h . Summarizing, the 
uncertainty in the measured metallicity of SNIa hosts is similar for 
the different metallicity tracers, e ~ 0. 1 dex, although it may be 
substantially larger for a given galaxy. 

We have drawn in Figs. [6] [7] and [8] an error bar indicative of 
the typical uncertainty in the measured host metallicity. These un- 
certainties are not large compared to the range of host metallicities 
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covered in the figures (horizontal axes). One can also compare the 
host metallicity uncertainty with the correction necessary to obtain 
Zi a , plotted in the same figures (vertical axes). Looking at the cor- 
rections deduced for the DTD of P08 (left panels of the same fig- 
ures) it can be seen that they are in general much larger than e. On 
the contrary, if the DTD of MlOb is used (right panels of Figs.[6]and 
0, AZi a is of the same order as e, especially for early-type galax- 
ies. Thus, if this last DTD were the correct choice, the error in the 
determination of SNIa metallicity would be dominated by the un- 
certainty in the measurement of the host metallicity rather than the 
difference between galactic and supernova metallicities. 



4 EXPLOITING GALACTIC DATABASES 

In this section, we explore the facilities provided by current 
galactic catalogues in order to obtain observational constraints to 
our theoretically derived MDF of SNIa. Ideally, one might use 
an observationally determined star formation history (SFH) of a 
galaxy together with a gas metallicity history, determined through 
metallicity-dependent synthetic stellar population (SSP) models, to 
recover a probability distribution function for the metallicity of 
SNIa progenitors. In practice, however, there are limitations due 
to both the discretisation of the catalogue in time bins of finite size, 
and the paucity of metallicities used in the SSP models. In Sec- 
tion !4. II we shortly address how to reformulate Eqs. |7]and|8]in or- 
der to calculate the probability that a SNIa progenitor was born in a 
given time bin. We apply then the new formalism to the VESPA cat- 
alogue, and discuss the results and the consequences of the paucity 
of metallicities in the SSP models in Section |4~2l We explain the 
main characteristics of the VESPA catalogue in the next paragrap h. 

VESPA ('VErsatile SPectral Analvsis'. lToieiro et alJkoOTo is 
a method to reconstruct star formation and metallicity hist ories 
from galactic spectra using SSP models. iToieiro et al.l d2009l) ap- 
plied their VESPA c ode to all galaxy spectra in the seventh data 
release of the SDSS dYork et aljEoOOl) . and compiled their results 
in a publicly accessible catalogue of stellar masses, star formation 
rates and metallicity histories of nearly 800,000 galaxies^. VESPA 
uses all of the available absorption features and the shape of the 
continuum (emission lines are not included in the analysis) in order 
to estimate the SFH of each galaxy, with variable time resolution, 
and the metallicity of the stars bom in each time bin, which traces 
the gas metallicity. The number of time bins depends on the qual- 
ity of the data on each galaxy. At the highest resolution, VESPA 
uses 16 age bins, logarithmically spaced between look-back times 
0.002 Gyr and 13.7 Gyr, however if the data are not of sufficient 
quality some bins may lack information. The VESPA catalogue 
contains results from several runs using different SSP models, IMF 
prescriptions, dust models, and galactic samples (for details see the 
above-mentioned papers from Tojeiro et al). The SSP models are 
supplied at metallicities of Z = 0.0004, 0.01, 0.02, and 0.04. 



4.1 Probability that a SNIa originated in a given time 
interval 

The present rate of SNIa explosions in a given host is the convolu- 
tion over time of the SFR and the DTD. The probability that a SNIa 
that explodes now (at time t) comes from a progenitor born at time 
f = t - r is given by Eq.|5] Assuming we know the DTD and that 



http://www-wfau.roe.ac.uk/vespa/index.html 



the SFR can be obtained from a galactic catalogue, both in discrete 
time bins, the probability that the SNIa was born at time bin i is (for 
a given galaxy): 



Pi 



(28) 



(29) 



where the addition in the denominator extends to all bins defined 
for that galaxy, and Mi is the mass of stars formed in time bin i. 

If the observational DTD is not known with enough resolu- 
tion (as compared to the time bin resolution in a given catalogue), 
one can turn to the model of P08 (Eq. |7J or its equivalent based 
on the DTD of MlOb (Eq.[SJ. What we actually need in order to 
use the galactic catalogue is the probability that a SNIa exploding 
now comes from a progenitor born within a time interval given by 
the temporal bins defined in the catalogue: say between look-back 
times Tj and r, + A?, : 

f T ' +A '' , (dM\ , 

P,<x J SQ-t'MM^—j dt', 

where we have used the DTD of P08. To use the catalogued data 
we transform the above integral using a mean SFR defined from the 
stellar mass formed in the time interval i and the duration of that 
time interval, Af,: 5, = M/Af,, and combine (dM/dt)df = dM: 

M f M(Ti) 

Pi oc — 1 <p(M)dM . (30) 

Af; Jikf(T;+Afj) 

The integral in the last equation can be easily computed given the 
IMF, the stellar lifetime law, and the temporal limits of each bin. 
Note that there is no need to know the proportionality constant for 
Pj nor even to re-normalize the IMF, because we know that J] Pi = 
1, if we add Pj for all the bins in a galaxy. Thus, the probability 
density can be normalized unambiguously. 

If the DTD of MlOb is adopted instead, Pj can be obtained 
from: 



Mj - 
Pj cc — -*P, 

Afj 



(31) 



where is the average DTD between t, and r, + Atj. 

From the previously determined probability that a SNIa comes 
from time bin i, Pj, and the catalogued metallicity for that same 
time bin, Z,, the metallicity probability density is given by: 



P'(AZ) = P'ilogZ, - logZ, IM ) = Pj . 



(32) 



This procedure provides us with a discrete set of AZ for which the 
probability density is known. 

4.2 Results using VESPA: dependence on the star formation 
rate and host age 

Here we show an example of the application of the formalism intro- 
duced in the last Section to real data. We have accessed the VESPA 
catalogue (runID=2) and classified the ga laxies as early or l ate- 
type according to the criterion proposed bv lDildavetailfcoirJ) on 
base of their u and r SDSS model magnitudes. Then, we have se- 
lected the galaxies with a minimum temporal resolution of three 
bins, in order to work with reasonable SFHs. Finally, we recovered 
90906 early-type galaxies and 40419 late-types. We stress that we 
have not filtered the catalogue in order to select which galaxies are 
potential SNIa hosts. From a theoretical point of view, every sin- 
gle galaxy can house binaries with the appropriate parameters (to- 
tal mass, secondary mass, initial separation, metallicity) to produce 
SNIa. Observationally it is another story, because surveys aiming to 
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detect SNIa have to pre-select target galaxies for which the obser- 
vations and measurements are more efficient. This difference has 
to been kept in mind when comparing between our theoretically 
derived SNIa statistics and the correlations found in observational 
campaigns. 

We have applied the formalism developed in the previous sec- 
tion (Eqs.l31landl32t to all the galaxies selected from VESPA in or- 
der to obtain the MDF of each galaxy. Unfortunately, at the present 
level of resolution of the catalogue the resulting MDF is of little 
practical use because of the limited spectra of metallicities of the 
SSP models. As the MDF is an statistical description of the dif- 
ference of metallicity between the host and the SNIa, the use of 
just four metallicities in the SSP models produces an unreliable 
extremely clumpy MDF. Keeping this in mind, we comment here 
briefly on the relationship between the metallicity correction, AZ, 
and the Zh ost obtained from the catalogue. Linear fits to these quan- 
tities are shown as dashed lines in Fig. [6] Given the uncertainties, 
it is striking the close coincidence between the points belonging 
to our theoretical models and the linear fits obtained with the cat- 
alogue, especially for late-type galaxies. The differences between 
theoretical models of early-type and late-type galaxies are more or 
less reproduced by the mean metallicities derived from VESPA. 

In view of the strong uncertainty affecting the metallicity 
statistics, we do not follow further with the analysis of the VESPA 
derived MDF, but we emphasize that the formalism presented in the 
previous subsection can (and should) be applied to future galactic 
catalogues based on improved data and metallicity resolution. In 
spite of the failure to derive useful MDF from VESPA, the data 
from the catalogue can still be used to show trends of the mean 
SNIa metallicity, Z 1;l calculated using the probability defined by 
Eq.[3T] vs. the galactic SFR and age. 

Figure [9] shows Zi a as a function of the sSFR (specific Star 
Formation Rate). We first notice the sharp cut on top of the points 
distribution at Z[ a = 0.4, which derives from the maximum metal- 
licity fed to the SSP models of VESPA. If we disregard this limita- 
tion of the models, we can still fit a tendency line to the data points, 
with the result that can be seen in the plot: SNIa exploding in high 
sSFR galaxies do have smaller Z [a than those that explode in passive 
galaxies. On the other hand, the tendency lines for early and late- 
type galaxies do not differ significantly. While the correlation is not 
strong, we note that a relationship between sSFR and SNIa lumi- 
nosity has been found in some obser vational studies of SNIa (e.g., 
to cite just one, I Sullivan et al .Il20l"ct who found that SNIa in low 
sSFR galaxies appear brighter on average than those in high sSFR 
galaxies, after applying stretch corrections). As a note of caution, 
we remark that our galaxy sample was drawn from the whole unfil- 
tered VESPA catalogue whereas each observational study of SNIa 
properties are based on particular selection function of hosts. Thus, 
a direct comparison with our results might not be meaningful. 

Figure [10] shows the distribution of Zi a as a function of the 
look-back time of the host galaxy (that is: f back = d/c, where d is 
the distance to the galaxy and c is the speed of light), that is itself a 
measure of the age of the universe at the local host time, and of the 
redshift. Besides the accumulation of points around Z = 0.04, the 
data show a tendency for smaller Z la at lower redshift. The trend 
of Zhost with look-back time (not shown in the Fig.) is opposite 
to that of Z Ia , i.e. older hosts are less metal-rich as expected from 
popular galaxy evolution models. The sign of the slope of Z Ia vs. 
fback is difficult to explain, as it implies that SNIa are on average 
less luminous at high reds hifts, contrary to w hat is derived from 
observational samples (e.g. iHowell et alj|2009l their Fig. 8). Any- 
way, we cannot derive any strong conclusion given the limitations 
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Figure 9. Mean Zj a as a function of the specific SFR, as derived from the 
VESPA catalogue (runID=2) and the DTD of MIOb. The solid lines are the 
results of linear fits to the points given by [Zi a ] = -0. 177 log (sSFR) - 2. 10 
(blue, early-type galaxies) and [Zi a ] = -0.196 log (sSFR) - 2.26 (red, late- 
type galaxies). 



of the present analysis, as the trend we detect might be a result of 
observational biases of the catalogued data. 

We have checked that the above trends are qualitatively inde- 
pendent of the assumptions concerning the DTD and the specific 
runs of VESPA that can be selected (galactic samples, SSP mod- 
els, IMF prescriptions, and dust models). We have also tested the 
scatter of the parameters of the tendency lines in Figs.l9landll0las 
a result of the uncertainties in the galactic properties stored in the 
catalogue (errors in the metallicity and star formation rates in each 
galactic bin). To this end, we have computed again the tendency 
lines for 200 random realizations of the galactic sets (early-types 
and late-types) by adding random noise according to the catalogued 
metallicity and SFR errors in every time bin of each galaxy of our 
VESPA sample. For each random realization, we have computed 
the tendency line and thereafter we have analyzed the distribution 
of slopes, their mean and standard deviation. The standard devia- 
tion of the slope of the tendency lines of sSFR and / back with respect 
to Zia is typically less than 0.002. 



5 CONCLUSIONS 

We have calculated the chemical evolution of homogeneous galaxy 
models together with the evolution of the supernova rates in or- 
der to evaluate the Metallicity Distribution Function of SNIa, 
MDF(AZ), i.e. the probability that the logarithm of the metallic- 
ity of a SNIa exploding now differs in less than AZ from that of its 
host. We have analysed several model galaxies aimed to represent 
from active to passive galaxies, including dwarf galaxies prone to 
experience supernova driven outflows. 

Our results show a remarkable degree of coincidence (in an 
statistical sense) between the mean Zi a and Zh 0S t- The dispersion 
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Figure 10. Mean Zi a as a function of the look-back time of the host galaxy, 
f back , derived from the VESPA catalogue (runID=2) and the DTD of Ml Ob. 
The solid lines are the results of linear fits to the points given by [Zi a ] = 
0. 147 f back -0.524 (blue, early-type galaxies) and [Z Ia ] = 0.064/ back -0.251 
(red, late-type galaxies). 

of Zi a in active galaxy models is quite small, meaning that Zh ost is 
a quite good estimator of the supernova metallicity, while passive 
galaxies present a larger dispersion. We have devised a procedure 
to correct the difference in metallicity between SNIa and their hosts 
(Sect. |3.2t , based on a linear fit given in Fig. [6] These results are 
insensitive to the choice of IMF, stellar lifetime and stellar yields, 
within the ranges we have explored. In contrast, the Delay Time 
Distribution of SNIa remains one of the main ingredients influenc- 
ing the difference of metallicity between SNIa and their hosts. 

We have discussed the use of different metallicity indicators 
(Fe vs. O, gas-phase metallicity vs. stellar metallicity). Using O (as- 
suming its evolution is representative of that of the CNO group) as 
a metallicity measure for late-type galaxies does not change appre- 
ciable the metallicity correction with respect to using Fe (assuming 
its evolution is representative of that of all metals It is remarkable 
that the metallicity correction for early- type galaxies (jar < 1 Gyr) 
is quite different for O than for Fe. On the other hand, the metallic- 
ities of SNIa are much better represented by those of the gas-phase 
than by the mean stellar metallicities. Thus, when possible, it is ad- 
visable to measure the host metallicity through that of its gas-phase. 

Finally, the results of the application of our formalism to a 
galactic catalogue (VESPA) suggest that SNIa come, in average, 
from smaller metallicity progenitors both at low redshifts and in 
galaxies with large star-formation activity. The paucity of metallic- 
ities used in the original stellar-population synthesis models used 
in the construction of the galactic catalogue does not allow us to 

5 The precise meaning of this statement is that the mass fraction of CNO 
elements in SNIa progenitors is represented by their mass fraction in the 
host galaxies, with similar accuracy as the whole metal mass fraction in the 
progenitors is represented by the whole metallicity of the hosts. 



build a significant MDF of SNIa based on the galactic histories 
stored in the database. However, in spite of the large uncertainty 
in the metallicity derived from the catalogue, the gross trends of 
Zi a vs. Zhost obtained from VESPA for different galaxy types are 
roughly consistent with our theoretical estimates. The derivation of 
improved, observationally based, MDF will be possible in the fu- 
ture if the SSP models include further refinements in their grid of 
metallicities. 

One of the main drawbacks of our approach is the use of ho- 
mogeneous one-zone galactic chemical evolution models. Actual 
galaxies are heterogeneous, and radial gradients of chemical com- 
position are routinely measured within them. A possible improve- 
ment of our model would b e to divide the galaxy in several non- 
interacting shells as in, e.g.. lMartinelli et ail ( ll998b . then comput- 
ing the galactic MDF accounting for the probability that a SNIa 
originates in each shell. Then, however, a complication would arise 
when comparing to observations of far away galaxies for which 
there is no possibility of measuring the metallicity gradient: what 
is the meaning of the galaxy metallicity in a context where each 
independent shell has its own chemical history? What is clear is 
that the derived MDF of SNIa would be characterized by a larger 
dispersion than the MDF we have obtained. In some sense, the true 
MDF of SNIa reflects a double dependence: in space, i.e. the posi- 
tion of the SNIa progenitor within its host galaxy, and in time, i.e. 
the moment in the galactic history when the progenitor was born. 
It is just this last dependence what we have studied in the present 
work. 
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